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Abstract 

This  research  project  was  initiated  to  explore  the  use  of  novel  analysis  techniques 
associated  with  emerging  nonlinear  methods  to  investigate  the  growth  of  capiUary-gravity 
waves  on  an  initially  calm  water  surface  following  the  abrupt  start  of  wind  forcing.  During 
the  course  of  the  research,  a  novel  wavelet-based  method  was  devised  and  implemented  for 
computing  the  local,  instantaneous  dispersion  relation.  The  method  is  completely 
quantitative,  and  conserves  energy.  While  averaged  results  showed  exponential  growth 
rates  that  agree  approximately  with  linear  theory,  it  is  clear  that  details  of  the  wave  crest 
distribution  that  are  important  in  determining  subsequent  nonlinear  developments  cannot  be 
captured  without  not  only  a  nonlinear  theory,  but  also  a  2D  one.  The  most  useful 
experiment  to  follow  up  this  work  would  not  be  the  obvious  increase  in  scale,  but  would 
raAer  concentrate  on  a  smaller  scale  facility  where  the  2D  wind  profile  could  be  very 
closely  controlled,  and  measured.  Only  then  can  the  forcing  effects  be  unambiguously 
separated  from  2D  nonlinear  growth  mechanisms. 
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Format  of  the  report 

Since  the  principal  findings  of  the  study  have  been  published  in  two  papers,  that  are 
attached  as  appendices,  the  report  itself  will  serve  mainly  to  put  these  in  overall 
context  of  the  stated  objectives,  and  to  make  some  general  conclusions  from  the  two 
years  work. 

Original  objectives 

Initial  wind  energy  input  appears  as  small  scale  capillary  gravity  waves,  whose 
evolution  and  transition  to  longer  wavelengths  is  thus  important  in  determining 
overall  air-sea  interaction  energy  budgets,  and  in  generating  surface  roughness 
signatures  that  are  detectable  in  SAR  imagery.  The  original  proposal  promised  to 
investigate  the  physical  mechanisms  of  this  process,  using  novel  experimental  and 
analytical  techniques,  combined  with  a  theoretical  analysis.  The  long-term  goal  is 
therefore  (a)  to  improve  basic  physical  understanding  of  small-scale  wave 
generation  mechanisms,  (b)  to  make  quantitative  predictions  of  their  role  in  ocean 
energy  fluxes,  and  (c)  to  apply  some  of  these  concepts  to  remote  sensing  design  and 
analysis. 

Experimental  and  Theoretical  Methods 
These  can  be  enumerated  as  follows: 

1.  The  Imaging  Slope  Gauge  (ISG)  measurement  technique,  together  with  the 
Complex  Wavelet  Transform  (CWT)  dispersion  relation  analysis  method  was  used 
to  generate  and  test  predictions  for  growth  rates  of  capillary-gravity  waves.  The 
ISG  method  was  adapted  and  refined  in  the  first  3-year  ONR-sponsored  project. 
The  CWT  method  was  an  entirely  new  method  that  came  from  the  same  effort.  An 
extensive  analysis  was  performed  on  data  sets  collected  at  the  Air-Sea  Interaction 
Research  Laboratory,  Wallops  Island.  Thorough,  quantitative  comparisons  with 
results  from  numerical  simulations  of  the  nonlinear  Schrddinger  (NLS)'  equation 
have  been  conducted. 

2.  Theoretical  models  incorporating  increasing  degrees  of  realism  ii.e.  ZD/ 
unsteady /nonlinear  effects)  were  developed,  following  indications  from  the 
results  of  #1  above.  In  particular,  the  theoretical  analysis  of  Caponi  ei  al  (1992)  was 
repeated,  and  a  quasi-steady  approximation  allowed  solutions  to  a  time  varying 
problem  to  be  estimated  by  integrating  over  a  boundary  layer  thickness  that  grows 
in  time  and  has  a  profile  according  to  experimental  measurements. 

3.  Certain  limits  were  identified  in  the  application  of  current  theoretical  models  and 
laboratory  experiments  in  dealing  with  real  problems  of  oceanographic  and  ONR 
relevance. 
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Brief  stunmary  of  results 

The  results  are  given  in  more  detail  in  the  attached  appendices,  but  the  most 

important  overall  features  can  be  summarised  here. 

1.  Experimental  results  The  initial  wavenumbers  and  growth  rates  of  the  most 
unstable  modes  are  well  predicted  by  a  linear  instability  of  coupled  water-air  shear 
flow  models.  The  initial  wavenumbers  are  very  close  to  theoretical  predictions, 
but  careful  analysis  of  the  growth  rates,  both  from  this  study,  and  those  in  the 
literature,  reveals  that  they  are  not  predicted  by  the  same  degree  of  certainty. 
Factors  of  2-4  difference  are  not  uncommon.  A  new  method  of  using  the  phase 
component  of  the  CWT  for  automated  recognition  and  tracking  of  defects  has  been 
developed.  This  'waveprint'  uniquely  characterises  a  wave  field  in  terms  of  its 
dislocation  sites,  or  critical  points,  for  a  given  wavenumber.  A  new  method  for 
estimation  of  friction  velocities  from  a  combined  CWT/ theoretical  model  has  been 
developed. 

2.  Theoretical  models  The  quasi-steady  analysis  shows  that  predicted  dispersion 
relations  are  very  sensitive  to  assumptions  about  the  growth  in  characteristic  scale 
of  the  air  and  water  boundary  layers.  A  simple  viscous  diffusion  model 
completely  fails  to  produce  order  of  magnitude  agreement  with  experiment.  If  the 
observed  boundary  layer  growth  is  specified  then  the  resulting  dispersion  relation 
differs  from  the  free-wave  problem  in  qualitatively  similar  ways  to  experimental 
data.  It  is  clear  that  a  more  complex  treatment  is  required.  Some  of  the  3D 
structure  in  the  surface  deformation  is  due  to  defect  patterns,  and  new,  analytical, 
defect  solutions  to  the  NLS  have  been  found. 

3.  Limits  The  current  state  of  the  theoretical  art  is  found  to  be  unsatisfactory. 
Existing  linearised  models  that  predict  growth  rates  and  mode  amplitudes  are 
very  sensitive  to  extra  inputs,  such  as  boundary  layer  thickness  and  its  time 
variation.  Absent  a  predictive  model  of  these  data,  even  for  controlled  laboratory 
conditions,  we  lack  a  predictive  understanding.  The  discrepancy  with 
experimental  data  becomes  more  acute  at  higher  wind  speeds.  The  problem  may 
indicate  the  inadequacy  of  the  2D  (or  ID),  linear  model  formulations.  None  of  the 
latter  stages  of  growth,  where  energy  moves  to  increasingly  lower  wavenumbers, 
are  satisfactorily  modeled  at  present. 

On  the  other  hand,  with  current  uncertainty  in  the  exact  form  of  the  wind  forcing, 
one  cannot  always  separate  the  contributions  from  non  uniformity  in  the  wind 
forcing  from  genuine  3D  instability  modes.  Paradoxically,  further  and  deeper 
understanding  of  wave  generation  in  ocean  conditions  may  depend  on  the 
development  of  new,  small-scale  laboratory  models,  where  the  precise  forcing  can 
be  known  and  varied. 


Implications 

Science/Physics:  The  current  theoretical  underpinnings  of  the  wind  wave  initialisation 
problem  have  been  shown  to  be  less  secure  than  previously  thought,  thanks  to 
quantitative  comparisons  of  unprecedented  scope  and  accuracy.  At  the  same  time, 
the  basic  mechanism  of  air/water  shear  layer  instability  for  initial  wave  generation  is 
very  solidly  supported,  and  future  theoretical  developments  must  be  based  on  this 
fundamental  mechanism.  The  existing  discrepancies  indicate  that  either  a  3D  and/ or 
a  nonlinear  model  will  be  required  for  further  progress. 
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Technology:  CWT  analysis  has  been  used  on  both  laboratory  and  ocean  data  now, 
since  the  instrumentation  of  the  latter  is  improving.  Technology  is  currently  being 
shared  between  the  two.  The  development  of  nonlinear  mechanisms  of  wave  growth 
occurs  rapidly  at  all  but  the  lowest  wind  speeds,  and  the  development  of  the  lower 
wavenumbers  that  persist  in  ocean  applications  is  not  yet  modeled.  The  new  friction 
velocity  estimate  has  been  developed  specifically  for  field  and  SAR  image  data 
applications. 

Collaborations  and  research  transitions  ^  xt  .  x 

The  current  project  is  a  collaboration  between  USC  (Spedding,  Browand,  Newton), 
NASA  WFF  (Huang,  Long)  and  Scripps  Institute  of  Oceanography  (Jahne,  Klinke). 
Internal  NASA  fimds  have  been  used  to  co-support  much  of  this  work.  An  AASERT 
studentship  applying  wavelet  transforms  to  oceanographic  measurements  was 
spawned  by  this  project  and  continues  at  USC.  The  collaboration  with  NASA 
Wallops  was  particularly  successful,  and  we  hope  for  further  co-operative  ventures. 

Since  the  WT  analysis  has  been  developed  and  conducted  on  the  dual  machines, 
purchased  with  NASA  funds  expressly  for  this  purpose,  the  technique  is  transitioned 
very  rapidly  to  the  NASA  government  domain,  and  therefore,  potentially,  to 
application  to  large  ONR-  and  NASA-sponsored  databases.  We  have  been  active  in 
ensuring  that  new  programs  are  documented  as  they  are  written  and  NASA  WFF  has 
a  complete,  documented  set  of  all  software  generated  under  this  agreement. 


Appendices 

A.  Spedding  GR  1996  2D  Complex  Wavelets  for  Analysis  of  Unsteady,  Wind-Generated 
Surface  Waves.  Nonlinear  Dynamics  of  Ocean  Waves.  Special  Issue  of  Naval  Research 
Reviews,  58(3) ,  30-39. 

B.  Spedding  GR,  Klinke  J  &  Long  SR  1996  Estimating  co(k)  in  an  unsteady,  wind-generated 
surface  wave  field  from  the  2D  complex  wavelet  transform  of  the  surface  slope.  In:  The  Air- 
Sea  Interface:  Radio  and  Acoustic  Sensing,  Turbulence  and  Wave  Dynamics.  Proc.  Int.  Symp. 
Marseille,  1993.  (eds  M.A.  Donelan,  W.H.  Hui  &  W.J.  Plant),  373-382. 
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Abstract 

Wavdets  and  their  transfonns  (WT)  have  numerous  and  diverse  applications,  and 
almost  as  many  shapes  and  sizes.  This  article  concentrates  solely  on  their  utility  as 
a  quantitative  meagre  of  spatially-localised  spectral  amplitudes  for  2D  data 
analysis.  The  {diysical  problem  concerns  the  time-dependent  growth  of  wind- 
driven  surface  waves,  starting  from  a  flat  surface.  Although  the  problem  lacks  the 
^  complexity  of  real  ocean  conditions,  it  is  argued  that  because  it  is  precisely 
denned,  a  combined  experimental  and  theoretical  program  can  advance  our 
I^ysical  understanding,  that  may  itself  be  projected  onto  field  data.  Certainly,  the 
advent  of  nonlineari^  in  the  wavefield  is  marked  by  intermittency  and  3D 
stmcture,  and  the  WT  is  an  ai^ropriate  basis  for  analysis.  The  wind-wave  problem 
will  be  used  as  an  example  where  the  simultaneous,  quantitative  measurement  of 
amphtude  and  i*ase  modulations  in  an  unsteady  wave  field  leads  to  estimates  of 
quantities  such  as  the  local  dispersion  relation  that  carmot  easily  be  measured 
otherwrse.  Furthermore,  the  lATT-derived  results  can  be  combined  with  simple 
analytical  models  to  measure  quantities  that  are  both  of  great  interest  to  the 
oceanographer,  and  that  have  traditionally  been  difficult  to  estimate.  Fmally, 
suggestions  are  made  for  future  work  and  take  a  surprising  direction. 


L  What  is  the  problem? 

Why  study  small-scale  wind  waves? 

For  the  purposes  of  this  article,  small-scale  waves  are  those  having  wavelengths  in  the 
cmtimeter  range.  Both  gravitatitmal  and  surface  torsion  effects  are  significant  and  so  the 
^ysical  proUem  is  an  interesting  balance  between  the  two.  From  a  practical  point  of  view, 
tltere  are  two  reasons  for  interest  Capillary-gravity  (eg)  waves  are  the  wave-like  disturbances 
that  first  grow  when  a  calm  water  surface  is  disturbed  by  wind.  Later  in  time,  the  energy- 
containing  scales  are  much  larger,  but  at  any  instant,  it  is  through  smaller-scale  motions  such 
as  eg  waves  that  the  instantaneous  surface  response  mediates  the  fluxes  of  momentum  and 
ener^  Mtween  the  atmosf^re  and  ocean.  TTie  second  practical  motivation  is  that  remote 
%nsmg  devices  such  as  SAR  are  sensitive  to  scattering  of  electromagnetic  radiation  primarily 
m  the  shorter  wavelengths  characterised  by  eg  waves.  It  is  only  by  modulation  of  the  shorter 
umvelengths  that  larger  scale  fe^res  are  observable,  (jenerally  it  is  sufficient  to  model  this 
effect  as  a  roughness’,  but  statistical  roughness  models  and  predictions  will  not  behave 
reahshc^y  unless  they  are  based  mi  sound  physical  principles.  As  we  shall  see  later,  the  eg 
wavefield  can  be  extrmnely  complex  in  structure,  and  further  understanding  leads  inevitably 
to  studres  that  fall  squarely  into  the  ‘Nonlinear  Dynamics’  theme  of  this  issue. 


Why  stu^  the  unsteady  problem? 

It  is  a  mild  understatement  to  observe  that  fluctuations  of  the  local  wind  stress  (in  both  space 
and  time)  ^r  over  the  ocean  surface.  While  steady-state  response  models  may  work 
^tisfactorily  in  some  statistical  sense,  the  transient  response  and  physical  understanding  of 
how  one  state  changes  to  another,  together  with  the  possible  interactions  between  a  time- 
y^mg  component  with  another  st^onary  one,  are  still  of  great  interest  Furthermore,  while 
the  luteansed  models  predict  the  initial  stages  of  the  surface  response  quite  well,  they  carmot 
predta  the  subsequettt  evolution  of  low  wavenumber  (long  wavelength)  modes,  3D  pattern 
tn^ihtres,  and  wave-wave  interactions  that  follow.  In  other  words,  there  still  is  a  prtrWem, 
and  rt  has  not  yet  been  solved. 


^ere  rs  a  second,  more  abstract  line  of  reasoning  that  also  supports  such  a  line  of  amuiry. 
^e  of  the  principal  achievements  of  nonlinear  dynamical  systems  theory  has  been  to 
dernonstrate  how  complex  {^nmnena  can  be  generated  by  simple,  deterministic  systems, 
wrth  few  degrees  of  freedom.  Hopes  have  been  raised  that  taslu  of  apparently  daunting 
complexity  may  yet  yield  to  analytical,  predictive  methods.  In  a  similar  vein,  starting  from 
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flic  veiy  evidoit  fart  that  the  real  sea  surface  state  is  extremely  complicated,  one  searches  for 
the  sirnfflest  possible  model  that  can  nevertheless  exhibit  sufficiently  rich  dynamics  that  it  can 
account  for  the  most  important  physical  mechanisms. 


Sta^g  ^m  scratch,  we  require  a  system  that:  (i)  is  clearly  defined;  (ii)  can  be  modeled 
analyhcally,  or  numerically;  (iii)  can  be  tested  experimentally,  and  (iv)  can  be  generalised  to 
situations  of  practical  interest  Condition  (iv)  is  rarely  satisfied  unless  (ii)  and  (iii)  involve  a 
degree  of  understanding  based  on  physical  principles,  and  the  failure  of  empiric^  matching 
of  model  parameters  and  coefficients  when  applied  outside  of  the  dmnain  for  which  they 
have  been  carefully  tuned  is  a  rather  familiar  problem  in  oceanography. 

The  spectnm  of  analyflcal  models  for  the  air-water  interface  is  bounded  on  the  one  hand  by 
sinall  amplitude,  two-dimensional  solutions  of  the  Rayleigh  equation  for  steady  mean  flows 
Math  boundaries  at  x  *  +/-«».  Correspondingly,  the  Blasius  solution  to  the  incompressible,  flat- 
plate  boundary-layer  equations  for  a  specified  origin  in  x  gives  the  fetch-limited  result 
Between  the  two  limiting  cases,  we  are  interested  in  the  time-dependent  problem,  in  the 
absoice  of  external  boundaries. 


One  might  questitm  the  relevance  of  laminar  model  flows  when  the  ocean  aimlications  are 
almost  always  turbulent  in  both  fluids.  However,  in  the  initial  transient  for  a  wind  profile 
staiwg  from  nwt  even  if  the  flow  itself  is  not  perfectly  laminar,  the  laminar  instability  modes 
comd  still  be  predominant  in  selecting  the  most  energetic  wave  modes.  In  this  respect  the 
laboratory  experiment  incorporating  some  degree  of  turbulence  in  the  growing  profiles, 
provides  a  bridge  between  model  and  reality. 

yfhy  bother  with  laboratory  experiments? 

In  the  previous  article,  Norden  Huang  has  considered  the  general  applicability  of  wave-tank 
ex^riments  to  real  ocean  conditions.  It  is  worth  reiterating  the  basic  philosophy  behind  our 
laboratory  work  and  its  relation  to  the  overall  problem  outlined  in  the  previous  paragraphs. 

In  contrast  to  fidd  da^  it  is  only  in  the  relatively  controlled  omditions  of  the  laboratory  that 
[ur^eters  can  be  v^ed  at  will,  independently,  and  repeatedly.  It  is  not  necessarily 
sigmficant  that  corditions  do  not  exactly  mat^  those  of  the  open  ocean,  because  the 
experiments  can  still  provide  non-trivial  tests  of  anal3^cal  or  conceptual  models  that  are 
themrelves  of  use  over  a  broader  range  of  parameters.  The  support  or  refutadmi  of  such 
mMels  can  still  provide  information  about  their  ability  to  predict  events  in  the  real  woild.  If 
and  when  the  model  fails,  it  can  be  extended  or  throMoi  away;  in  either  event,  the  conditions 
under  which  it  fails  are  rxiw  knowru 

In  the  remainder  of  this  article,  a  closely-linked  theoretical  and  experimental  research 
program  will  be  d^nbed,  where  wavelet  analysis  will  be  used  as  a  tool  for  identifying  current 
hmitatKNis  and  pointing  to  future  requirements. 


2^  The  contribution  of  wavelet-bascd  methnri^ 

Wavelet  as  atomic  units 

At  the  time  of  writing,  ^re  are  11  texts  on  wavelets  in  our  university  bookstore,  and  some  are 
8®*^'  For  a  concise  and  readaUe  introduction,  it  is  hard  to  beat  the  Algorithms  and 
Applications  book  by  Meyer .  A  more  comprehensive  and  mathematical  overview  is 
presented  in  the  intr^uctory  volume  of  the  series  on  wavelets  edited  by  Chui*,  and 
geophysical  Allocations  of  varying  quality  can  be  found  in  a  later  volume*  in  the  same  series. 
Farge  has  wntten  a  fine  review  on  wavelets  and  turbulence,  and  further  collections  of 
theoretical  advaiKes  and  aj^licatitms  are  accumulafing.*** 

^e  of  the  great  strengths  of  the  field  is  in  the  mathematical  depth  and  rigour  of  much  of  the 
denv^on  and  analysis,  which  has  revealed  connections  with  hitherto  unrelated  techniques, 
from  Litflewood-Paley  analysis  of  the  1930’s  to  modem  quadrature  mirror  signal  processing 
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algoritluns  and  spline  functional  analysis.  Nevertheless,  a  zero-th  order  degree  of 
sophistication  in  understanding  wavelet-based  signal  analysis  will  suffice  for  our  puq^res. 

On  this  level,  wavelets  are  simple  oscillatory  functions  that  have  compact  support,  aiw  mus 
can  be  used  for  local  space-scale  decompositions.  The  original  function,  fiX)  is 

expanded  as  a  function  of  an  additional  scale  variable,  a 

fix)  fia.x).  (1) 


The  decomposition  is  effected  by  a  convolution  of  fix)  with  multiple  copies  of  a  single 
wavelet  function,  y(jc),  that  are  systematically  shifted  and  dilated  to  cover  a  range  of 
characteristic  scales,  a,  and  positions,  b. 

Provided  certain  admissibility  criteria  are  met  (y(x)  must  have  finite  energy,  zero 
mean...consult  refs.(l,2,4)),  then  the  choice  of  y(Jt)  is  somewhat  ^itrary.  It  is  possible  to 
maife  the  strong  assertion  that  the  physical  worid  is  much  more  likely  to  be  comi»sed  of 
collections  of  wavelet-type,  compactly  supported  functions,  rather  thm,  say,  infiratmy- 
extended  oscillatory  functions  as  in  Fourier  series.  It  follows  that  it  is  only  reasoMble  to 
decompose  real  world  observations  (and  computations)  onto  such  atomic  rmits.  Farge  ,  for 
example,  draws  attrition  to  the  fact  that  a  decomposition  of  a  complex  fluid  flow  into  its 
elementary  coherent  structures  (if  such  an  object  exists),  is  better  approximated  by  a 
functional  base  with  compact  support  than  the  usual  Fourier  base.  Alternatively,  one  might 
search  for  basis  functions  whose  properties  are  extracted  from  the  data,  or  from  kr»i^ 
solutions  to  the  particular  governing  equations.  There  will  likely  be  future  progress  in  this 
direction,  but  it  is  in  fact  unnecessary  to  adhere  to  strictly  to  the  atomic  construct. 


Wavelets  as  tools  for  quantitative  description 

The  keyword,  of  course,  is  ‘construct’,  and  one  is  free  to  constract  a  wavelet,  or  any  other, 
basis  to  represent  the  worid  in  many  different  ways.  For  example,  one  of  the  pnncipal 
advantages  of  Fourier  representations  is  the  considerable  mathematical  convenience  m 
computation  of  PDF’s,  having  little  or  nothing  to  do  with  any  implied  physical  ^is. 

Similarly,  provided  one  is  aware  of  the  properties  of  the  decompositiori  and  its  role  in 
determining  (even  distorting)  the  interpretation  of  the  world,  one  can  still  use  a  wavelet  basis 
for  descriptive  purposes.  This  is  not  all  the  same  thing  as  stating  that  wavelet-atoms  actually 
exist,  or  that  they  have  any  special  dynamical  significance.  Hence,  we  simply  state  the 
following; 

1.  The  objective  is  to  make  quantitative  measurements  and  interpretations  of  local  energy 
densities  in  a  2D,  time-varying  surface  wave  field. 

2.  The  most  straightforward  way  to  do  this  is  through  continuous,  complex-valued  wavelets. 

The  definition  of  the  wavelet  transform  (WT),  described  in  words  above,  can  be  stated  more 
precisely: 

fia^b)  =  J ^ 

The  transfonn  is  a  sum  of  the  convolution  of  the  translates  and  dilates  of  y  (x)  withy(x)  (the 
*  denotes  the  complex  conjugate).  In  one  dimension,  a  useful  choice  of  y(x)  is  the  Monet 
wavelet 

=  (3) 

As  Figure  1  demonstrates,  it  is  a  complex-valued  function,  and  the  7c/2  phase  shift  tetween 
real  and  imaginary  parts  allows  the  modulus  of  its  coefficients  to  be  fiw  from  oscillations  due 
to  changes  in  phase  in  an  otherwise  constant  amplitude  signal.  The  wavelet  indicates 
the  {^ase  in  the  instantaneous  frequency,  or  wavenumber,  associated  with  a  particular  scale,  a. 
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ko  is  constant  (t3rpically  e  ]S.S,6{)  so  there  is  a  fixed  number  of  cycles  in  an  analysing 
wavelet,  regardless  of  a.  Consequently,  the  spatial  resolution  improves  with  increasing 
wavenumber,  or  decreasing  a.  llie  way  in  wMch  it  does  so  varies  automatically  with  a,  and 
the  balance  between  resolution  in  wavenumber  space  and  in  physical  space  is  in  some  sense 
an  op^um  for  the  Morlet  wavelet,  within  the  overall  constraints  of  Heisenberg’s  uncertainty 
principle. 

The  WT  conserves  energy,  and  a  form  of  the  Parseval  relation  can  be  written, 

—  ® 

where  is  a  normalisation  constant  that  depends  mi  the  form  of  y(x).  The  right  hand  side  is 
an  integi^  of  a  function  that  is  proportional  to  the  q;>atial  distiibutimi  of  energy  density  per 
unit  sdde.  It  is  possiUe  to  integrate  contributions  fr^  all  x  to  arrive  at  a  globid.power 
spectrum. 


£(<J)  =  ;r  (5) 

E{fl)  is  a  Gaussian-smoothed  version  of  a  power  spectrum  that  could  equally  well  have  been 
obtained  from  Fourier  transforms.  Of  course,  in  case  we  are  not  oblig^  to  perform  the 
summation  over  x,  and  the  point  is  to  note  that  the  energy  densities  are  precise,  quantitative 
measures  of  localised  spectral  components  (provided  that  there  ate  no  artifacts  from  the 
mutual  interaction  of  overlap^ng,  non-orthogonal  wavelet  coefficients). 

The  extension  to  two  dimensions  of  this  basic  approach  is  fairly  straightforward  and  details 
can  be  found  in  ref.  9.  As  before,  we  require  a  complex  wavelet  function,  and  for  real-valued 
signals  (2D  surface  slope  measurements,  for  example)  it  is  sufficient  to  cover  the  hidf-plane  in 
wavenumber  space.  Now  that  y=y(x,y),  the  wavelet  can  be  designed  so  that  it  can  be  either 
directional-specific,  or  non-directional.  A  non-directional  wavelet  responds  equally  to 
wavevectors  at  all  angular  orientations,  and  in  two  dimensions  such  a  Action  covers  an  arc 
that  sweeps  througli  180"  in  [k„  k,)  (Fig.  2). 

Arc  and  its  close  relatives  form  the  core  of  the  data  analysis,  giving  a  measure  of 
instantaneous,  wavenumber  energy  densities  and  wave  p^ase  as  a  ftmction  of  {x,y}. 


3.  An  integrated  anaivtical/experimcntal  program 

Theoiy 

The  time-varying  solution  for  perturbations  to  the  interface  between  two  viscous  fluids  with 
specified  shear-flow  (xofiles  can  be  solved  numerically  by  integrating  coupled  Orr- 
Sommerfeld  equations  for  air  and  water.  Both  Valenzuela**  and  Kawai“  have  done  that,  two 
decades  or  so  ago.  Exponential  growth  rates  for  eg  waves  were  predicted  (and  measured,  in 
the  case  of  Kawai)  and  compared  reasonably  well  with  microwave  backscatter  measurements 
of  Larstm  &  Wri^L“  This  body  of  worlc  represents,  still,  the  state-of-the-art  in  the  numerical 
analysis  of  the  linearised  problem. 

The  stability  of  piecewise  linear  shear  profiles  in  air  and  water  was  recently  examined  by 
Caponi  et  al.^  One  can  take  their  steady-state  model  and  integrate  forwards  in  time  by 
allowing  a  boundary  layer  thickness  to  grow,  thus  generating  a  quasi-steady  solution.  The 
sensitivity  of  the  solution  to  the  assumed  growth  of  the  boundary  layer  thickness  will  tell  us 
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ofiuch  about  the  behaviour  of  the  numerical  solutions,  and  the  directions  for  further 
improvements. 

Experiment 

Experiment^y,  we  investigate  the  growth  of  disturbances  on  an  initiaUy  calm  water  surface 
when  the  wind  starts  from  rest,  and  at  sufficient  fetch  so  that  upstream  boundary  conditions 
can  be  ignored  (at  least  in  the  water).  The  configuration  is  a  reasonably  close  experimental 
approximation  of  the  theoretical  approaches  discussed  above. 

The  measurements  are  inade  with  a  modified  form  of  the  Imaging  Slope  Gauge  (originally 
described  by  JMhne  &  Riemer*^,  which  uses  surface  refraction  of  intensity-graduated  incident 
.  light  m  derive-ibe  instaritaiieous  along- ^  cross-wind  slope  over  the  field  of  view.  s,(x,y,r)  or 
Syix.yfy  can  be  reranstructed  according  to  the  directitm  of  the  light  intensity  gradient  The 
surfa(%  slope,  is  directly  related  to  the  wave  energy,  and  the  WT  time  series,  'sJix,yA,t) 
contains  information  about  the  local  energy  densities,  and  about  the  local  wavenumbdr- 
frequency  relation.  Model  predictions  can  thus  be  carefully  checked.  For  the  first  time,  the 
unsteady,  local  dispersion  relation  can  be  measured,  and  related  to  the  2D,  instantaneous 
structure  of  the  wavefield. 

Matching  the  two 

Neither  existing  models,  nor  practicable  experiments  can  provide  complete  and  generalisable 
information  about  the  time-varying,  3D,  w&vt  and  velocity  fields.  It  will  emerge  that  ttie 
theoretical  models  are  very  sensitive  to,  and  experiments  mostly  ignorant  of,  the  detailed 
shape  and  time-evolution  of  the  water  and  air  boundary  layers.  However,  in  combinatkm, 
mutual  consistency  checks  and  judicious  matching  of  empirical  data  with  analytical  relations 
can  lead  to  new  methods  for  measuring  and/or  inferring  the  values  of  physical  variables 
(surface  stress,  boundary  layer  height),  not  only  from  laboratory  data,  but  also  from  remote 
sensing  measurements.  Two  examples  of  this  kind  of  reasoning  will  be  given.  In  both  cases, 
the  glue  that  holds  the  pieces  together  is  the  Arc  wavelet  transform. 


^  The  time  evolution  of  capillarv-gravitv  waves 

The  surface  wave  growth  for  a  mean  wind,  f/,o  =  6.8  m/s  (measured  at  10  cm  above  the  calm 
water  surface)  is  shown  in  Hgure  3.  The  time  series  proceeds  from  left  to  right,  top  to 
totmrn,  and  is  for  wind  forcing  from  the  right  The  dimensionless  time  x  =  UrotPC,  where  X  is 
the  fetch,  ranges  from  3  to  5.5,  so  the  upstream  origin  is  always  felt  in  the  air  flow,  even  at  the 
early  times.  The  group  velocity  of  the  initial  water  waves,  however,  is  much  smaller.  At  50 
cm/s,  for  example,  24  seconds  are  required  to  travel  the  fetch  distance  of  12m. 

It  is  clear  that,  even  at  early  times,  the  observed  wavefield  is  never  perfectly  two-dimensional. 
The  grey-scale  encoding  of  the  slope  magnitude  gives  tiie  initial  disturbarKes  an  organ-pipe 
^^arance,  but  they  always  have  Y-shaped  ‘defects’,  or  ‘dislocations’,  where  successive 
waves  become  entangled.  In  fact,  it  was  the  presence  of  the  defect  patterns  that  first  drew  our 
attention  to  the  wind-wave  problem,  as  we  wondered  vriiether  there  may  be  simUarities  with 
other  defect  structures  in  fluid  turbulence.’*’” 

Even  though  the  waves  are  not  perfectly  two-dimensional,  calculation  of  the  global  power 
spectral  densities  and  their  rate  of  change  with  time  shows  that  the  most  amplified 
wavenumbers,  and  their  initial  exponential  growth  rates,  can  be  quite  well  predicted  by  linear 
theory  (Figure  4).  That  is  ^e  for  the  first  five  images  (the  top  row)  in  Figure  3.  During  the 
cour%  of  about  1.25s,  the  initially  most-amplified  disturbances  have  reached  their  maximum 
^plitude,  and  in  the  top  right  panel,  the  waves  look  as  two-dimensional  as  they  ever  will  be. 
Shortly  afterwards,  the  surface  patterns  become  much  more  complicated  as  the  initial  modes 
saturate  arul  further  growth  occurs  in  sidebands  of  bofli  higher  aixl  lower  wavenumber.  By 
the  end  of  the  series  longer  wavelength  disturbances  predominate.  A  convincing  physical 
explanation  for  this  well-known  phenomenon  has  never  been  given.  Note,  however,  that  it 
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does  ^ipear  to  be  a  duration-limited  event,  since  information  propagated  by  water  waves  from 
upstream  boundaries  will  not  yet  have  reached  the  measurement  window.  Further  insights  can 
be  derived  from  wavelet  analysis. 


5.  Wavelet.based  analysis 

Local  energy  density  distributions 

The  remailcs  concerning  the  global  evolution  of  instability  modes  are  based  on  Fourier 
analysis.  The  2D  Fourier  transform  can  provide  a  perfectly  good  description  of  the  average 
luoperties  of  the  surface  dope,  and  has  very  good  spectral  resolution.  Having  identified  the 
wavenumbers  of  interest,  it  is  now  possible  to  use  the  localised  properties  of  tlw  wavelet 
transform  to  relate  the  s^tral  densities  to  the  structure  of  the  wavefidd. 

The  ^tid  distribution  of  growth  of  the  initially  most  amplified  mode  at  {/,o  =  6.8  m/s  can  be 
seat  in  Figure  5.  The  ener^  density  distributions  are  quite  patchy,  and  remain  so.  By 
analogy  with  defects  in  mixing  layers,  where  the  defects  are  the  sites  for  the  first  (and 
subsequent)  subharmonic  transitions,  it  is  reasonaUe  to  ask  whether  the  globally-measurable 
sideband  instabilities  ate  located  at,  or  near,  defect  sites  in  the  surface  wave  field.  The  open 
circles  denote  foe  location  of  defects,  identified  by  eye  in  the  equivalent  time  series  of  surface 
slq;)e.  Again,  by  ^ogy  with  foe  mixing  layer  results,  the  local  energy  density  at  the 
fiirklamental  scale,  Oq,  would  be  a  minimum  here,  while  the  sideband  modes  would  peak  at  foe 
same  location.  At  first  glance,  the  circles  do  seem  to  be  consistemly  beside,  rather  than  on, 
the  highest  energy  densities  in  Oo.  However,  Figure  6,  which  is  the  same  plot  for  foe  lower 
wavenumber  sideband,  a.i,  shows  that  the  circles  are  not  associated  with  local  maxima  at  this 
scale.  It  can  also  be  found  that  the  other  sideband,  characterised  by  wavelet  scales  a„  is 
correlated  quite  well  with  maxima  in  u.,,  and  so,  also  not  with  Oo.  The  high  arxl  low 
wavenumter  modes  correlate  well  with  each  other,  but  wily  weakly  with  minima  in  the  initial 
tnode.  Initially,  the  defects  do  not  appear  to  [day  a  significant  rdle  in  the  development  of 
sidebands.  Although  orie  cannot  expect  ^atial  correlations  to  be  maintained  in  a  dispersive 
wavefield,  foe  variation  in  phase  speed  with  wavenumber  is  rsfoer  weak  at  this  particular  range 
of  scales,  and  the  wave  pattern  defects  or  dislocations  do  not  appear  to  be  involved  in  the 
growth  of  sidebands. 

Pattern  dtfects  and  the  waveprint 

The  correlation  of  pattern  defects  with  various  postulated  transitions  arid  instabilities  could  be 
made  more  precise  if  a  reliable  method  for  identifying  the  defect  sites  were  available,  arui  we 
recently  found  one. 

Recall  that  it  is  ^cifically  defects  at  a  particular  scale  that  we  would  like  to  investigate,  aiKi 
the  WT  can  be  employed  as  a  kind  of  filter  so  that  only  contributions  from  this  scale  are 
counted.  Images  sudt  as  Figures  5  &  6  are  clearly  of  i»  use,  but  the  phase  of  the  WT  is. 
Hgure  7  is  a  Mruuised  image  of  foe  phase  of  the  WT  at  scale  Not  only  can  one  confirm 
the  defect  locations  that  were  originally  identified  by  eye  (circles),  but  it  is  also  possiUe  to 
follow  their  time  evolution.  Ntuie  of  ^  defect  tt^logies  in  foe  circles  remains  fixed  over 
time;  instetKl  connections  with  jneceding  or  following  waves  are  initiated  and  brcAen.  In  a 
similar  spirit  to  foe  descriptitm  of  fluid  flows  by  critical  poirtts,"  the  wavefield  topology  can 
be  character!^  by  foe  type  and  location  of  its  defects,  together  with  the  distribution  of  the 
highest  amplitude  wavelet  coefficients.  Just  as  fingerprints  can  be  categorised  by  their 
defects,  so  the  wavefield  can  be  described  by  its  waveprint  One  can  look  for  particular  defect 
structures  in  an  ap{»opriate  analytical  model,  and  ^auf  Newton  at  USC  has  made  a  start  in 
describing  a  new  famUy  of  defect  solutions  to  the  norilinear  Schrddinger  equation.'* 
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Estimation  of  local  dispersion  relation 

The  WT  ^ves  a  measure  of  the  spatially-localised  energy  densities  as  a  function  of 
wavraumber,  and  s^e  the  me^rcments  are  also  well-resolved  in  time,  the  local  frequencies 
can  be  calculate  from  ©  =  d^/dt.  This  is  possible  because  the  complex  WT  also  givM  the 
instant^ous  phase  distribution,  ^(x,y),  associated  with  each  scale  (as  shown  in  the 
wavepimt).  Figure  8  straws  how  the  calculation  proceeds.  The  flowfield,  u(x,y),  that  maps 

9\X.y4i)  to  energetic  scales,  a  (associated  with  physical  wavenumber,  k),  is  used 

to  compute  the  local  frequency  as 

to(j:,y)  =  |A:(jr,y)jM(jr,y).  (6^ 

TOere  is  M  comparable  time-de^ndent,  2D  theory  that  predicts  such  a  detailed  dispersion 
relation,  but  values  over  aU  {x.y)  can  be  collected  together  and  compared.  Figure  9^  an 
ex^ple,  where^the  solid  curves  have  been  drawn  over  an  arbitrary  range  of  k  from  the  result 
of  C^m  et  al  for  the  mtnnsic  ©(A:)  of  small  perturbations  on  a  surfa^  between  two 
piecewise  luie^  shear  profiles.  The  experimental  results  differ  from  the  theoretical  curves, 
wrach  IS  Mt  alro^ther  surprising  since  the  wavelengths  (and  frequencies)  that  actually  have 
TC  Ivgest  amplification  rate  in  practice  will  depend  on  the  time  history  of  the  growth  in 

tetwm  the  <o(t)  ajKl  the  solid  cuive.  either  the  surface  drift  or  the  boundary  layer 

^  mferred.  This  is  an  example  of  a  combined  experimental/theoretical 
modelmg  for  simultaneous  measurement  of  a  variable  whose  estimation  could  be  of  great 
importance  m  the  field.  * 

Estimation  of  surface  stress  from  lab.  or  field  data 

let  us  sui^rase  tl^  is  krrawn  (from  existing  experimental  data  or  field  measurements. 
The  maxim^  dope  ^ckness,  D,  can  be  written  as  tiie  ratio  of  the  mean  wind,  (7,  to  the 
gradient  of  Widy  at  the  water  surface  .  uki 


fnetion  yelwity, ,  is  a  measure  of  the  wind  stress  exerted  at  the  water  surface,  which 
also  vanes  with  aUldy: 


air  Vist^ity  and  density,  respectively.  Combining  the  two  expressions, 
the  friction  velocity  can  be  related  to  D,  ©  k  , 


2  vU 

“■‘T' 


vi^ity.  An  esttnute  of  D  is  therefore  equivalent  to  an  estimate 
of  the  fn^on  vdocity.  To  see  how  this  relationship  can  te  incoiporated  into  the  overall 

scheme  is  the  combined  theoretical  model  for 
^*)  together  with  expenmental  measurements  of  the  same  quantity.  Given  a  measured 
dispersion  relation,  and  a  selection  criterion  for  choosing  one  particular  curve,  such  as  the 


8 


maximum  ampUfication  rate,  then  the  difference  between  measured  and  predicted  curves 
can  be  attributed  to  the  surface  drift  current,  The  corrected  ratio  of  the  w^r  and  air 
velocities  are  returned  to  the  theoretical  model,  and  upon  convergence,  a  predict^  'I 
D=D  can  be  used  to  measure  u..  It  is  noteworthy  that  the  measurements  do  not  have  to  come 
fromllaboratory  data.  In  principle,  a  time  series  of  SAR  images  could  provide  the  same 
information,  because  the  WT  can  give  estimates  of  <o(A:)  from  these  data  also. 

The  two  examples  above  demonstrate  how  the  newly-acquired  ability  to  me^ie 
instantaneous  mnpUtude  and  phase  information  for  local  spectral  com^nents  c^cad  to 
novel  and  potentially  powerful  measurement  techmques  that  have  application  m  ootn 
laboratory  and  ocean  data. 


7.  Concluding  remarks 

The  wind-wave  startup  problem  has  been  explored  as  an  exainple  of  the  application  of 
wavelet-based  measurem^t  and  analysis  techniques  to  uns^y.  inhom^eneous 
the  extension  to  these  domains  that  brings  us  closer  to  realishc  ^an  randibons. 
into  nonlinear  regimes  that  have  just  begun  to  be  investigated.  The  advanced  measurement 
techniques  have  application  to  botii  laboratory  and  field  conditions. 

With  reeaid  to  the  particular  wind-wave  problem,  it  is  obvious  that  the  wavefield  is  ^ways 

quite  complex  in  structure,  but  the  dynamical  significant  of  tlj  X 

^ious.  In  attempting  to  distinguish  between  structure  that  is  due  to  .  , 

the  water  surface  deformation,  or  that  merely  reflects  imprecision  m  sp^fication  jjf 

forcing,  we  need  a  smaller-scale,  more  easUy  controll^le  expenmenL  The 

wind  orofile,  and  its  time  history,  must  be  specified.  Then,  an  experiment 

modtJations  in  both  along-  and  cross-wind  directions  to  gener^ 

would  a  wind-forced  equivalent  of  existing  systematic  investigations  of  liquid  surface 

instabilities.**^ 

Paradoxically,  the  conclusion  that  future  experiments  ought  to  be  in  smaller,  nrorc  precisely- 
conliolled  conuineis.  moves  sUU  liiither  from  tnie-Wue  "“““SVh?.  wtem  Rey^s 
numbers  are  large,  and  many  simultaneous  processes  com^te  an^or  mteract  wth  e^h  otne  . 
Recalling  the  introduction  however,  it  is  likely  dial  the  basic  instab^ttes  md  ^ 

mechanics  wiU  only  be  isolated  and  identified  with  a  more  selective,  and  P^>sely  specified 
experiment  We  hope  to  report  on  the  results  of  similar  experiments  m  the  future. 
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Figure  captions 

Figure  1  Is  the  world  composed  of  wavepackets  like  this?  The  Morlet  wavelet  is  a 
smooth,  continuous  wave  function,  modulated  by  a  Gaussian  envelope.  Real  and  imaginary 
parts  are  shown  in  continuous  and  dotted  lines.  Whether  or  not  it  is  a  tnie  atomic  unit  of  the 
{^ysical  world,  it  can  still  be  a  convenient  function  basis  for  general  description. 

Figure  2  The  2D  wavelet  function  Arc  is  defined  most  conveniently  in  Fourier  space 
where  for  a  given  scale,  a,  it  can  be  seen  to  interrogate  an  arc  th^  covers  all  wavevectors  of 
magnitude  lU,  regardless  of  their  orientation. 

Figure  3  The  time  evolution  of  alongwind  surface  slope  following  a  change  in  mean 
wind,  f/,0,  from  0  m/s  to  6.8  m/s.  The  change  in  wind  ^ed  is  aj^roximately  linear  over  5 
seconds.  The  time  series  shown  here  begins  at  ts5.3s,  i.e.  after  Uxo  bas  reached  its  steady 
value.  The  field  of  view  in  each  image  is  approximately  14  x  11.5  cm. 

Figure  4  Experimentally-measured  initial  exponoitial  growth  rate,  p,  vs.  wavenumber, 
(symbols  connected  by  solid  lines),  for  f/,o  from  2.8  -  10.2  m/s.  The  solid  symbols  are 
from  experiments  of  Kawai,  and  the  dashed  lines  are  results  from  Kawai’s  Orr-Sommerfeld 
model  integration  for  parameters  that  closely  matched  those  ex^riments  ((/,o  is  given  in 
brackets).  At  moderate  windspeeds,  the  agreement  is  good,  but  is  very  sensitive  to 
assumptions  about  how  the  boundary  layer  grows,  and  so  remains  somewhat  empirical, 
because  there  is  no  analytical  model  for  the  boundary  layer  growth.  A  simple  difhision 
model  does  not  work  at  all,  and  a  pessimist  mig^t  wdl  argue  that  the  agreement  simply 
signifies  that  both  parties  have  done  similar  experiments.  As  selected,  the  best  match  occurs 
at  moderate  winds^eds,  but  the  extremes  are  not  so  well  covered. 

Figure  5  The  modulus  of  the  wavelet  transform,  IWTI^,  at  the  initially  most  amplified 
scide  (ko=4.0  rad/cm,  cm).  The  colour  bar  mailing  is  constant  over  the  series,  which 
covers  the  exponential  growth,  saturation,  and  decline  in  energy  densities  around  this  scale. 
Although  the  wavenumber  and  growth  rates  are  quite  well  predicted,  as  described  in  Figure  4, 
the  distribution  is  far  from  uniform  in  {jc.y}. 

Figure  6  IWTI..,  for  same  time  series  as  Figure  5.  Maximum  growth  of  the  lower 
wavenumber  mode  (k.,=2.2  rad/cm,  X.,=2.9  cm)  occurs  around  the  middle  of  the  third  row, 
just  when  the  peak  of  the  initially  most  amplified  mode  in  Fi^re  S  begins  to  saturate  and 
decay.  The  boxes  marked  ‘k,’  are  located  at  the  energy  density  peaks  and  correlate  well  with 
peaks  at  the  higher  wavenumber  mode  (k,=7.9  rad/cm,  Xo=0.8  cm;  not  shown). 

Figure  7  The  waveprint  The  sign  of  the  phase  of  the  WT  at  allows  defect  sites  to  be 
located  with  little  ambiguity,  enabling  automated  identification  arul  tracking  algorithms  to 
woric  reliably. 

Figure  8  How  to  calculate  the  instantaneous,  local  dispersion  relation  o)(A:,  x,  y)  from  a 
time  series  of  complex  wavelet  transforms.  Given  the  surface  slope  distribution,  s/x,y)  (top 
row),  at  two  closely-spaced  times,  r„  tj,  the  wavelet  transform  of  each  can  be  separated  into 
modulus  and  i^uise  components  (rows  2  &  3).  For  all  scales  having  significant  energy 
densities,  the  flow  field  that  maps  the  isophase  lines  from  their  position  at  ti  to  their  position  at 
t2  can  be  computed  from  an  optical  flow  technique.  Note  that  frequencies  higher  than  the 
apparent  Nyquist  limit  can  be  found  because  the  structure  in  the  phase  field  ^ows  large 
displacements  to  be  computed  in  a  global  optimisation  solution.  Finally,  only  measurements 
that  occur  where  IWn  >  IWTl*^  are  retain^. 

Figure  9  An  example  comparison  of  measured  atik)  with  theoretical  predictions  in 
solid  curves.  Each  symbol  marks  the  average  value  of  the  collection  of  points  at  a  particular 
wavenumber  where  IWTI  exceeds  a  threshold  value.  Each  spatial  location  contributes  a  dot. 
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^ H  line  whose  length  shows  the  variation  in  measured 

^  '*  adequately  described  by 

iine»  theory.  After  t=4.95,  the  scatter  increases  and  by  t=5.23,  there  are  no  wavelet 
coefficients  that  exceed  the  threshold  value  at  the  originally  most-amplified  wavenumber. 

^  procedure  for  estimating  wind  friction  velocity  from 
Tit  ^  ®‘dier  laboratory  or  field  data. 

quantities  in  water  and  air,  respectively,  and  subscripts  m  and  d 
Identify  measured  and  predicted  values.  manup 
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Abstract.  The  initial  formation  and  growth  of  capillary-gravity  waves  under  wind-induced  sur¬ 
face  stress  is  studied  by  measurement  of  the  instantaneous,  2D  distribution  of  the  along  wind  slope 
at  a  fixed  fetch.  A  new  technique,  using  information  from  the  complex  2D  wavelet  transform, 
is  described  for  the  calculation  of  the  local,  instantaneous  dispersion  relation.  The  technique 
is  validated  on  simple  test  fields.  The  dispersion  relation  measured  from  the  wave  tank  data 
is  consistent  with  predictions  from  linear  stability  theory  at  early  times,  but  departures  can  be 
noted  later  on.  Some  physical  interpretations  are  offered,  and  possible  nonlinear  mechanisms  are 
discussed. 

Key  words:  wavelet  transform,  dispersion  relation,  nonlinear  waves 


1  Introduction 

Classical  linear  stability  analysis  has  been  partially  successful  in  predicting  the 
generation  and  growth  of  the  initial  wave  field  on  a  calm  water  surface  following 
an  impulsive  start  in  the  wind  stress.  However,  neither  the  Phillips  resonant 
interaction  [1,  2],  nor  the  Miles  critical  layer  instability  mechanism  [3,  4]  has  been 
able  to  incontrovertibly  account  for  laboratory  and  field  measurements  of  growth 
rates  and  wave  initialisation.  Kawai  [5]  made  careful  measurements  of  the  initial 
growth  of  wind  generated  surface  waves,  and  combined  this  with  a  stability  analysis 
by  numerical  integration  of  the  Orr-Sommerfeld  equation  with  initial  and  boundary 
values  suggested  by  the  experimental  data.  It  was  found  that  linear,  shear  flow 
instability  mechanisms  could  successfully  predict  growth  rates,  phase  speeds,  and 
wavenumber-frequency  relations  in  the  observed  wave  field,  to  within  experimental 
accuracy.  It  was  further  pointed  out  that  the  subsequent  evolution  of  long-crested 
to  short-crested  waves  could  not  be  expected  to  be  governed  by  linear  mechanisms, 
as  indeed  was  indicated  by  ’complicated  movements  of  spectral  peaks’.  These 
spectra  were  quasi-steady  estimates,  computed  by  the  maximum  entropy  method 
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(c/[6])  on  data  segmented  in  time. 

Following  recent  advances  in  local  space-scale  decompositions  by  means  of 
wavelet  functions,  it  is  now  possible  to  perform  a  systematic  quantitative  anal¬ 
ysis  of  such  unsteady  phenomena,  without  some  of  the  limitations  imposed  by 
the  Fourier  decomposition,  and  where  the  balance  between  resolution  in  physical 
and  spectral  space,  scales  automatically  with  wavenumber,  within  the  limits  pre¬ 
scribed  by  the  uncertainty  principle.  The  objective  of  this  paper  is  to  describe 
the  development  of  a  method  for  the  measurement  of  spati^ly-localised  frequency 
and  wavenumber  relations  of  unsteady  wind- generated  surface  waves,  and  to  give 
examples  of  its  initial  application  to  recent,  2D^  instantaneous  wave  slope  data. 
This  will  include  the  later  times  when  nonlinear  mechanisms  may  be  anticipated 
to  be  significant. 


2  Experimental  data 

The  experiments  were  conducted  at  the  NASA  Wallops  wind-wave  interaction 
facility,  which  is  described  in  detail  in  [7].  The  aspects  of  significance  to  this  paper 
are  shown  in  Fig.  1.  Air  is  drawn  by  suction  over  a  75  cm  deep  water  channel 
over  a  working  section  of  approximately  18m  length  and  0.9m  width.  The  air  gap 
over  the  still  water  surface  is  46  cm.  Observations  of  the  stream  wise  surface  slope 
were  made  at  a  fetch  of  6.5m,  using  an  Imaging  Slope  Gauge  (ISG)  with  a  field  of 
view  of  approximately  18.5  x  14cm.  The  single  experimental  run  treated  here  was 
made  with  a  mean  reference  velocity,  11$  —  5.2  m.s“^  (measured  5cm  above  the 
surface,  and  15  cm  upstream  of  the  ISG).  Digitising  of  the  video-encoded  CCD 
signal  began  7s  after  the  fan  was  turned  on,  at  which  time  U$  =  5.0  m.s”^.  By 
frame#50,  8.6s  after  fan-on,  the  steady  value  of  11$  had  been  reached.  The  a;,  y  and 
axes  of  a  Cartesian  coordinate  system  run  streamwise,  spanwise,  and  vertically 
up  from  the  water  surface,  respectively. 

The  physical  problem  is  thus  somewhat  complicated  not  only  does  it  involve 
temporal  and  spatial  growth  of  viscous  boundary  layers  in  both  water  and  air, 
but  the  initial  conditions  and  start  up  are  much  less  homogeneous  in  space  and 
instantaneous  in  time,  than  one  would  perhaps  prefer,  and  careful  physical  inter¬ 
pretations  must  be  made  in  this  light. 

The  ISG  works  on  principles  described  in  [8].  A  strong,  homogeneous  (in  y) 
light  source  illuminates  a  rectangular  box  filled  with  particles  in  suspension  that 
scatter  the  light  so  that  its  intensity  decreases  with  downstream  distance.  In  the 
absence  of  any  surface  deformation,  light  rays  pass  normal  through  the  surface  and 
the  CCD  illumination  is  also  homogeneous  in  j/,  monotonically  decreasing  in  x.  A 
nonzero  surface  slope  in  x  {s^)  causes  light  rays  to  be  refracted  and  arrive  at  a 
different  x  location  on  the  sensor  array.  Thus,  given  the  intensity  distribution  on 
the  array,  it  is  possible  to  compute  the  deflection  and,  hence,  the  value  of  at  each 
position  in  {x^y}.  With  careful  calibration  procedures  [9],  this  calculation  can  be 
performed  so  as  to  reduce  errors  due  to  nonlinear  effects  to  less  than  0.05|5a:l. 


paper2.tex  -  Date:  October  31,  1994  Time:  14:59 


Dispersion  relation  in  unsteady  wind  waves 


3 


Fig.  2  shows  a  time  sequence  of  over  a  1.5  second  period,  during  which  time 
the  mean  and  maximum  slopes  increase  by  more  than  a  factor  of  10.  Even  at 
early  times,  the  wave  field  is  never  strictly  one-dimensional,  or  uniform  in  y.  It  is 
unclear  as  to  how  much  this  reflects  the  inhomogeneity  of  the  original  forcing.  At 
frame#70,  y)  reaches  its  most  ordered  state,  in  which  the  spanwise  correlation 
length  is  a  maximum.  After  this,  the  spajiwise  correlation  length  decreases,  and 
shortly  after  the  end  of  this  sequence  the  wave  field  becomes  quite  disordered,  as 
noted  originally  by  Kawai  [5],  and  more  recently  in  [10]. 


3  2D  Wavelet  functions  for  2D  data 

A  very  readable  introduction  [11]  and  review  of  wavelet  applications  [12]  have 
recently  become  available,  and  the  specific  design  and  application  of  2D  wavelet 
functions  for  analysis  of  2D  real  signals  is  discussed  in  [13].  Denoting  /w(a,6),  as 
the  wavelet  transform  (WT)  of  the  function  /(x),  a  is  the  analysing  wavelet  scale, 
and  resonance  at  any  particular  value  of  a  can  be  related  to  a  wavenumber  k  in 
the  signal.  This  information  is  localised  at  a  shift  from  the  origin,  6,  equivalent  to 
physical  location  x. 

It  is  convenient  to  think  of  the  wavelet  transform  as  a  product  in  Fourier  space 
of  the  wavelet  function  g{ak)  with  the  Fourier-transformed  signal,  /(^),  and  [13] 
introduced  the  Arc  wavelet  as  a  smooth  function  for  general  use  in  a  complex, 
non-directional  specific,  energy-conserving  WT  for  2D  data  analysis.  Here,  the 
alongwind  surface  slope  is  =  dsfdx^  and  since  Sx{x,y)  also  contains  contribu¬ 
tions  from  Sy{x^y)  projected  on  a:,  the  response  ought  to  be  weighted  by  cos0, 
where  6  is  the  angle  between  k  and  the  mean  wind,  i/co-  We  therefore  introduce 
the  Cos  wavelet  function: 


k  e  :  g{ak) 

k  E  :  g{ak) 


g(ak)  =  cosSexp 

0, 


(1) 


where  0  =  {kx>,k).  It  describes  a  cosine-weighted  arc  in  wavenumber  space  with 
Gaussian  cross-section,  centred  at,  and  shifted  from  the  origin  by  \ko\/a,  with 
width  l/(cra).  Isometric  surfaces  of  the  Cos  function  in  wavenumber  and  physical 
space  are  shown  in  Fig.  3.  Owing  to  the  smooth  decay  in  5cos?  9Cos  has  improved 
localisation  in  physical  space  over  the  original  gArc^ 

The  Cos  WT  can  be  computed  at  each  timestep  from  y),  and  the  result  can 
be  expressed  as  the  amplitude  |/w(fl,6)|,  and  the  phase,  Z/w(a,6)  of  the  complex 
WT. 
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4  How  to  calculate  u{k)  from  fw{a,b,t) 

Since  one  has  local  phase  information  from  the  complex  WT,  then  the  vector  field 
u  that  performs  the  mapping  of, 


(2) 

where  /.fw(a,b)  is  known  at  two  closely- spaced  timesteps,  can  be  used  to 

estimate  a  local  frequency.  Smce  u  is  measured  at  fixed  wavelet  scale,  a,  it  is 
associated  with  wavenumber  |^1,  and 


w  =  |A|u. 


(3) 


Various  methods  can  be  devised  for  computing  u.  Although  details  will  not  be 
given  here,  a  hybrid  scheme  using  a  novel  combination  of  optical  flow  and  cross¬ 
correlation  methods  was  found  to  work  reliably.  [14]  can  be  consulted  for  overviews 
of  each  of  these  classes  of  methods  separately.  The  2D  spatial  cross- correlation 
seeks  to  find  the  local  displacement,  r,  at  discrete  interrogation  windows,  such 
that 

H  f{x,h)f{x  -1-  r,t2)  (4) 

x^R 


is  maximised,  and  then 


u  =  fl(t2-ti). 


However,  simple  solutions  to  this  probl^  blow  up  in  the  presence  of  the  /oca/ phase 
discontinuities  of  ±x  found  in  Z/vy(o,6),  and  it  is  desirable  to  find  some  globally 
smooth  solution,  such  as  commonly  computed  using  optical  flow  techniques. 

u  was  computed  using  a  modification  of  the  technique  of  [15]  who  kindly  made 
available  their  original  codes.  We  proceed  by  minimising  the  cost  function. 


e(“)  =  +  7ies(u)  +  72  (V.M)^}</fi. 


(5) 


e/(u)  and  es{u)  are  terms  related  to  the  L2  error,  and  the  smoothness  (both 
global  terms),  respectively,  and  71  and  72  are  weightings  on  the  smoothness,  and 
a  divergence- free  constraint  that  was  introduced  by  [15]  for  fluid  flow  applications. 

When  the  displacement  field  due  to  u  is  large  compared  to  the  grid  resolution, 
a  multiresolution  relaxation  method  is  used  to  avoid  convergence  on  local  mini¬ 
ma  in  eq.(5).  Beginning  from  the  original  data,  successive  binary  expansions  are 
performed  until  the  local  displacements  do  not  exceed  one  mesh  size.  Given  an 
initial  estimate  of  the  flow  field,  the  solution  at  any  level  n  is  interpolated  out  to 
level  n  -fl  and  used  as  the  initial  approximation  for  the  n  +  I*'*  iteration.  Circular 
Gaussian  functions  were  used  for  the  pyramid  expansion,  and  a  conjugate-gradient 
method  was  used  for  the  iterative  solution  of  successive  linear  approximations  to 
eq.(5). 
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The  choice  of  weights  71 ,2 >  was  made  empirically  following  extensive  numerical 
experiments  on  synthetic  wave  fields.  These  test  also  showed  that  the  multireso¬ 
lution  optical  fiow  (MROF)  method  could  not  be  relied  upon  to  give  the  correct 
absolute  value  for  a  constant,  smooth  displacement  field,  over  a  large  range  of 
displacements,  for  fixed  71,2*  This  is  problematic,  since  the  variation  of  u  with  |A:1 
is  of  particular  interest.  A  convenient  and  reliable  solution  turns  out  to  be  to  use 
a  cross-correlation  method  in  order  to  generate  the  initial  first  estimate  for  the 
lowest  resolution  optical  flow  computation.  Owing  to  the  grid  expansion,  the  local 
behaviour  of  eq.(4)  is  not  a  liability,  and  one  immediately  benefits  from  a  very 
accurate  first  estimate.  The  optical  flow  computation  then  preserves  an  expan¬ 
sion  of  this  initial  estimate  while  maintaining  a  71-weighted  global  smoothness. 
Further  numerical  experiments  confirmed  the  robustness  of  this  hybrid  technique 
(XMROF)  over  a  range  of  displacement  magnitudes, 

A  block  diagram  summarising  the  calculation  methods,  from  initial  digitised 
data  to  a;(fc)  is  given  in  Fig.  4.  The  local  frequency,  u  is  computed  from  eq.(3) 
for  all  |/w(a,6)|  >  |/w(a,  The  threshold  value  was  set  to  80%  of  the 

maximum  value  of  |/w(fl,?)|  at  each  timestep. 

Extensive  numerical  simulations  were  performed  on  synthetic  data  fields,  and  on 
Lfw{ctjb)  derived  from  real  data,  but  given  a  known  form  of  u{k).  These  results 
will  be  reported  in  detail  elsewhere.  In  general,  the  mean  difference  between 
spatially- averaged  estimates  of  dJ(fc),  and  the  true  value,  wa>s  always  less  than  5% 
in  u  and  less  than  10%  in  \k\.  In  conditions  closely  matching  the  experimental 
signal  to  noise  ratios  and  resolution  in  x  and  the  mean  errors  were  about  half 
that. 

6  Time  dependent,  local,  dispersion  relation  for  wind-generated 

surface  waves 

Although  measurements  are  available  for  u;(|fc|,x,y),  they  will  be  averaged  over 
{x,y}  and  plotted  only  for  |a;|  for  brevity  and  ease  of  comparison.  This  dispersion 
relation  is  plotted  for  3  timesteps  in  Fig.  5(a-c).  Fifteen  frames  separate  each 
plot,  which  therefore  represents  data  from  a  completely  new  wavefield  that  has 
convected  into  the  measuring  location.  The  measured  data  points  can  be  compared 
with  the  dispersion  relation  that  can  be  derived  from  the  linearised  momentum 
and  continuity  equations  at  the  interface,  as  given  by  Caponi  et  a/.[16], 

w  =  y/i9k{l  -  a)  +  ak^),  (6) 

where  a  is  the  ratio  of  air  to  water  density  (=  1.21  x  10”^),  g  is  the  acceleration 
due  to  gravity,  and  <7  is  the  surface  tension  (74  cm^s“^).  The  measured  phase 
speeds  are  doppler-shifted  by  the  water  drift  velocity  at  the  surface,  Uo-  Since 
we  lack  accurate  simultaneous  measurements  of  Uq,  a  constant  offset  is  added  to 
place  the  solid  curve  on  the  data  point  occuring  at  the  predicted  most  amplified 
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wave,  at  ko  =  3.63  rad.cm"^  (An  independent  measurement  of  the  wavenumber  of 
the  most  amplified  wave  from  the  spatially- averaged  centroid  of  the  IVTT’p  gives 
A:q  =  3.67  ±  0.05  rad.cm"^).  This  offset  gives  a  drift  current  of  14cm. s”^,  which 
is  consistent  with  previous  results  reported  for  this  facility  [17],  but  at  .027?7oo  is 
significantly  below  the  .05  or  .04i7oo  sometimes  assumed  in  the  literature. 

The  measured  a;(|fc|)  in  Fig.  5a  agrees  with  the  theoretical  prediction  to  within 
experimental  uncertainty  over  the  whole  range  of  |fc|  where  there  is  significant 
energy.  At  the  next  timestep  (Fig.  5b),  the  range  of  energetic  |fc|  has  contracted 
around  ko,  consistent  with  qualitative  observations  made  earlier  for  Fig.  2  at 
t=65/70.  The  agreement  with  the  linear  theory  over  this  range  remains  good.  At 
t=:80,  the  range  of  energetic  modes  is  significantly  expanded,  with  energy  having 
appeared  in  both  lower  and  higher  wavenumbers.  It  is  greatly  reduced  around 
ko-  Above  ko,  eq.(6)  predicts  the  data  well,  but  there  is  increased  scatter  and 
difference  between  the  two  at,  and  below  ko. 

6  Discussion  &  Conclusions 

The  simple  relationship  of  eq.(6)  from  UD  linear  theory  appears  to  hold  well  until 
the  latter  stages  of  the  sequence  dicussed  here,  despite  the  fact  that  the  wavefield 
is  never  strictly  UD{i.e.  spanwise  uniform).  Nevertheless  the  initial,  most  unsta¬ 
ble  wavelength  is  well  predicted.  Caponi  et  at  [16]  include  a  discussion  of  the 
time- varying  interpretation  of  their  analytical  results,  which  were  obtained  from 
piecewise  linear  approximations  to  the  wind  and  water  velocity  profiles,  although 
the  solution  itself  did  not  involve  time- dependent  growth  of  these  profiles.  Quite 
specific  predictions  are  made  concerning  the  growth  rates  and  most  unstable  modes 
at  different  wind  speeds  and  we  look  forward  to  making  detailed  comparisons  over 
the  same  parameter  range.  The  results  of  Kawai  [5]  thus  far  appear  to  be  confirmed 
in  the  current  study,  and  will  also  be  compared  further. 

At  the  later  stages  corresponding  to  the  presumed  onset  of  nonlinear  effects,  we 
note  here  growth  of  modes  of  both  higher  and  lower  wavenumber  than  ko.  This 
is  illustrated  in  Fig.  6,  where  it  appears  that  growth  of  these  modes  occurs  in 
closely  associated  spatial  locations.  Although  an  earlier,  preliminary  wavelet  anal¬ 
ysis  [10]  indicated  that  significant  energy  developed  in  subharmonic  wavelengths, 
the  improved  resolution  and  quantitative  accuracy  of  the  current  data  allow  the 
distinction  to  be  made  that  the  majority  of  the  sideband  growth  actually  occurs 
in  bands  centred  at  |fco  ^ko-  The  sum  of  these  is  2ko,  consistent  with  the  res¬ 
onance  frequency  condition  for  a  Benjamin- Feir  type  instability.  However,  Bliven 
et  al.  [18]  suggested  that  this  mechanism  is  suppressed  by  active  wind  forcing, 
and  it  is  unlikely  to  be  significant  at  this  scale  range.  There  is  an  alternative,  sim¬ 
ple  physical  interpretation,  that  asymmetric  steepening  and  growth  of  single  wave 
crests  is  manifested  cts  different  modes  which  would  then  necessarily  be  correlated 
in  physical  space. 

In  this  case,  the  natural  tendency  to  think  in  terms  of  Fourier  modes,  or  their 
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spatially  localised,  wavelet-based  cousins,  may  be  unnecessarily  restrictive.  It 
could  prove  interesting  instead  to  allow  the  wavelet  function  basis  itself  to  auto- 
adapt  with  time,  similar  in  spirit  to  the  work  of  Coifman  and  colleagues  [19,  20]. 
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7  Figure  Captions 

Figure  1.  Simplified  schematic  of  the  wind- wave  generation  facility.  The  x’s  mark 
approximate  locations  of  pitot  tubes,  and  a  capacitance  wire  probe  for  wave 
elevation  is  marked  by  an  ’o’.  The  fan  control,  ID  digitising  of  pitot  tube  and 
wave  height  data,  and  sampling  of  the  2D  video  signal  are  synchronised  using 
3  personal  computers  and  one  Sun  workstation. 

Figure  2.  Time  sequence  of  surface  slope  values  computed  from  the  ISG  mea¬ 
surements,  numbered  in  frames.  Since  the  sampling  rate,  determined  by  the 
video  framing  rate,  is  30Hz,  the  total  sequence,  occupying  45  frames,  covers 
1.5  seconds.  The  doppler-shifted  group  velocity  is  such  that  an  entirely  new 
wave  field  is  covered  every  15  frames.  The  color  bar  is  scaled  to  cover  the 
maximum  and  minimum  slope  symmetrically  about  zero  at  the  last  time  step 
in  the  sequence. 

Figure  3.  The  Cos  wavelet  function  in  Fourier  space  (a),  and  its  real  and  imagi¬ 
nary  parts  (b)  &  (c)  in  physical  space. 

Figure  4.  Block  diagram  summary  of  the  method  for  calculating  d;(|it|)  from 
Currently,  the  minimum  At  is  determined  by  the  video  framing 
rate.  The  two  video  fields  are  treated  identically  on  either  side  of  the  line 
of  symmetry  in  the  block  diagram,  until  the  XMROF  method  combines  both 
sides  to  calculate  IWT^  i->  LWTt^.  All  boxes  marked  with  an  asterix  include 
state-of-the-art  advances  in  analytical  techniques  from  the  last  2  years. 

Figure  5.  Evolution  of  the  dispersion  relation  for  the  wind-generated  surface  wave 
field  of  Fig.  2.  (a)  t=50,  (b)  t=65,  (c)  t=80.  Diamonds  are  the  mean  measured 
values  at  each  |fc|.  The  solid  line  is  from  eq.(6),  drawn  over  the  whole  range 
of  k  for  which  the  WT  was  computed. 

Figure  6.  Isosurfa,ce  of  |/H^(a,6)|,  thresholded  at  50%  of  the  maximum,  for  t=80. 
Surface  intensity  variations  are  only  to  assist  in  showing  the  3D  shape.  The 

vertical  scale,  a  runs  from  Xmax  =  4.6cm  at  the  bottom,  to  A„,i„  =  0.6cm  at 
the  top. 
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